Bound on quantum computation time: Quantum error correction in a critical 

environment 
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We obtain an upper bound on the time available for quantum computation for a given quantum 
computer and decohering environment with quantum error correction implemented. First, we derive 
an explicit quantum evolution operator for the logical qubits and show that it has the same form 
as that for the physical qubits but with a reduced coupling strength to the environment. Using this 
evolution operator, we find the trace distance between the real and ideal states of the logical qubits 
in two cases. For a super-Ohmic bath, the trace distance saturates, while for Ohmic or sub-Ohmic 
baths, there is a finite time before the trace distance exceeds a value set by the user. 



Introduction- All physical implementations of quantum 
information processing face the inexorable reality of be- 
ing embedded in an environment that causes decoher- 
ence [l|. There are many strategies for dealing with this 
fact 0, quantum error correction (QEC) being one of 
the most versatile [11 . As QEC will likely be used in 
any large-scale implementation of quantum information 
processing, it is important to define and to quantify its 
limits. 

There have been several recent advances in under- 
standing the limits of fault tolerant quantum computing 
[3] • Part of this effort has been on proving "resilience" , 
the notion that any desired accuracy of quantum com- 
putation may be attained by concatenating levels of the 
QEC code y\. In particular, arguments for resilience in 
correlated environments have been constructed either by 
using techniques based on sums over faulty paths d, [7| 
or by reducing the problem to an almost stochastic one 
through scaling Q . In this Rapid Communication, we fo- 
cus on a related question: Given a certain quantum com- 
puter and an environment, what is an upper bound on the 
time available for computation? We provide an answer 
for a broad class of environments using a Hamiltonian 
formulation, including those where correlation effects are 
induced by gapless modes (i.e. critical environments). 

The main results of our argument are as follows. First, 
the coarse-grained quantum evolution of logical qubits 
is essentially the same as that of physical qubits, up 
to a renormalized coupling to the environment. Thus 
QEC maps our generic environment-computer interac- 
tion model onto itself, a property which has proven useful 
for other ways of suppressing decoherence such as dynam- 
ical decoupling |9|]. Second, we use this coarse-grained 
quantum evolution to find the maximum time available 
for computation, as assessed by the trace distance be- 
tween the real and ideal states of the computer. There 
is a regime where computation can continue indefinitely 
and so is resilient, while in other regimes the maximum 
time depends strongly on the QEC code, number of log- 
ical qubits, and environment-computer interaction. 

Environment- computer interaction and hypotheses- 



Many physical noise cases can be directly modeled by 
the "spin-boson model" [lo| . However, as originally 
proposed pT| . this model has a much more fundamental 
appeal. 

Let us assume that the qubits are already under the 
protection of an initial strategy, such as a decoherence 
free subspace or dynamical decoupling Though it 

deals with the dominant effect, it is unlikely to solve the 
decoherence problem completely. The inevitable residual 
interaction between the computer and the environment 
imposes a pointer basis for the qubits, which we call the 
z direction of each individual qubit. Another reasonable 
assumption is that the environment consists of a very 
large set of quantum degrees of freedom with some in- 
trinsic dissipative mechanism. Hence, in the absence of 
the qubits, the environment will be in a local minimum 
of its energy landscape. Our next (crucial) assumption is 
that linear response describes the influence of the com- 
puter on the environment. In that case, the computer- 
environment interaction cannot take the environment out 
of its local minimum, and so we may use the harmonic 
approximation to describe the quantum fluctuations of 
the environmental degrees of freedom. By construction, 
then, the environment-computer interaction term is lin- 
ear. Thus, we use the well-known and experimentally 
relevant spin-boson model to discuss the limits on pro- 
tection that QEC and fault-tolerant methods can yield. 

Having thus settled the model of the strongest channel 
of residual decoherence, we consider the presence of an 
additional weaker transverse channel, denoted by x. If 
the first channel were absent, the preceding argument 
applied to the transverse channel leads to a model of 
the same form but with a much weaker coupling: Aa, <C 
Xz- Hence, we describe the residual decoherence of the 
quantum computer by two bosonic baths (h — 1), 



Ho 



= {x,z} |k|^0 



(1) 



The fla k obey standard commutation relations, and 
LOa.k = i-^a (|k|/fco)^° , where ko and ujq are constants with 
dimensions of momentum and frequency, respectively. 
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and dynamical exponent. The environment- 

computer interaction has the form 



Hi= ^/"w^ <^ 



(2) 



a={x,z} X 



where cr^'^ 
sition X, :: 



are the Pauh matrices for the qubit at po- 
stands for normal ordering, and : /"(x) : = 



a.k 



H.c 



with 



(|k| /fco) which defines the exponent Sq,. Here, the 
environment has spatial dimension D and smallest mo- 
mentum 27r/L, and is a constant with dimensions of 
momentum. There is no a priori restriction on including a 
third bath; however, it would be a redundant description 
of the possible errors. All the bosonic averages performed 
are done with respect to the bosonic vacuum with no ini- 
tial entanglement between the computer and the bath. If 
there were some initial entanglement, it could be mod- 
eled using a finite temperature in the bosonic correlators, 
thus introducing an exponential scale. Since our goal is 
to calculate an upper bound for the computational time, 
we do not consider this case. 

In conjunction with this model, we make some assump- 
tions about the computer and the QEC method, (i) 
Gates are fiawless and are done much faster than the 
QEC period A. (ii) State preparation and measurements 
are done perfectly, (iii) Lowest order perturbation the- 
ory in Hi is enough to describe the evolution during a 
QEC step, (iv) All the syndromes indicate a non-error 
result; that is, we consider the most favorable quantum 
computer evolution, as any other evolution will involve a 
larger leak of information to the environment . 
Uncorrectable errors and the quantum evolution- The 
first step of the quantum calculation is to consider the 
evolution operator in the interaction picture [//(A, 0) — 

Tt exp —iJ^dtHj(t) up to a time A ~, just before error 

correction is applied. QEC divides errors into classes that 
can be distinguished from each other; however, within 
each class, different errors are not differentiated by the 
syndrome. For each logical qubit the syndrome breaks 
the evolution into m(A+,0) = X^i^o^'*^*' where i indexes 
the N possible syndromes of that qubit Each one 
of these evolution operators has a "good" and a "bad" 
part: = + f^i'^'' ^ where represents a logi- 

cal error. These logical errors are uncorrectable (or "bad 
evolutions" p|). 

Following our hypotheses, within the QEC period A, 
we approximate the evolution operator by expanding to 
lowest order in the couplings Aq. (Technically, the ex- 
pansion parameters are AqA and not simply Aq.) Thus, 
for a single qubit, 

Uj{A,0)^l-z ^A,A:r(x,0):aS. (3) 

a={x,z} X 

For a code of distance dc, one finds that the lowest order 
term that must be kept in each logical qubit is of order 



dc in the coupling to the environment. For illustration, 
consider the smallest distance-3 code, namely, the 5-qubit 
code 0. At the end of a QEC period, there are 4^ pos- 
sible configurations for the five qubits. They are divided 
into 4^ groups with distinct syndromes; however, each 
group has 4^ elements that cannot be distinguished by 
the code. We choose to analyze the evolution for which 
all the syndromes are the "no error" type. This yields 
the quantum evolution operator 



a/3 A A 



7;o(A,0) « l + zA3 Y 

x.,ct,p ,k 

X :r(x„0)::/^(x„0)::/'^(xfe,0):a^. 



(4) 



with X labeling the logical qubits, a,l3 — {x,z}, and 
j, J, k — {1, 5} labeling the physical qubits inside the 
logical qubit x. Each coefficient rj^-^ has two possible 
values, 77f|4 = ?7f|5 = ?7ff4 = rifi^ = ri^l^ = r]ll^ = rjlf^ = 
'^Ifs = '7523 = "nlib — 1 ^i^d zero for all other indices. 

Now, we use the commutation relations of the free 
bosons to normal order the evolution operator in Eq. Q , 

z;o(A,0)«i + * ^/,ffcA„A:r(x„0): 

X [ap,k + (XpA)' :/^(x„0)/^(xfe,0):] a^, (5) 

where aaij = (AqA)^ X^k^^o 1^^",^!^ "^xp [-ik • (x^ - x^)]. 
Equation ([5]) is written for the 5-qubit code with no 
concatenation; it is straightforward to generalize it to 
a larger distance or concatenated code. For instance, 
the level-1 concatenated code requires 25 physical qubits 
with the coefficients rj changing accordingly. In this case, 
A includes the time needed to extract all syndromes (in- 
cluding level-1 syndromes), and uncorrectable errors ap- 
pear at higher order in Aq. 

The evolution operator for a logical qubit can be 
abridged by rewriting Eq. ([5|) as 



wo(A,0) Kl + iA 



E 



(A;+r„):r(x,0): a«, (6) 



where x is the average position of the physical qubits be- 
longing to the logical qubit. A* = Aa X)/3,i,j,fc 
is the effective coupling constant, and rQ(x, 0) = 
AaE/5,.,,,fc'7r/fc(A/3A)2 ://'(x„0)/''(xfc,0): accounts for 
higher-order corrections. The latter dresses the single 
logical qubit amplitude a^jk by the interactions with 
other logical qubits. If the spatial separation of logi- 
cal qubits is at least S while that of the physical qubits 
within a logical qubit is ^, then Fq generates corrections 
of order (^/S)'''^" to observable quantities, where Sa is 
the smallest scaling dimension of the /„. For simplicity, 
we assume that ^ ^ S; hence, since we are seeking an 
upper bound on the computing time, we can ignore the 
Fq corrections. 

Another scenario to consider is when the physical 
qubits do not interact with each other, ^ — > cx). In this 
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case, apjk —>■ and the leading correction will come from 
Fa ■ This demands a slightly different organization of the 
argument: It leads to a different definition of the effective 
coupling constant but does not imply that there are no 
"uncorrectable errors" (see, e.g., Ref. [13!| for the stochas- 
tic error model). Most of the following discussion can be 
readily adapted to this case following arguments similar 
to those in Ref. tSjj which we therefore omit here. 

The steps outlined earlier result in the following quan- 
tum evolution operator for the logical qubits: 



(7) 



As a direct consequence of the coarse graining used in 
Eq. ([3]), note that the ultraviolet frequency cutoff is A^^. 

Equation ([7]) is a remarkable expression: It shows that 
in the long wavelength limit the logical gubits obey the 
same dynamics as the physical qubits. In other words, 
QEC maps the "spin-boson" decoherence model onto it- 
self. There are, of course, several ways to reduce A* : (i) 
engineer the position of the physical qubits, (ii) change 
the distance of the code, or (iii) concatenate the code. 
Nevertheless, as long as A* and Ta are not strictly zero, 
there will be degradation of the information in the logical 
qubits. Thus, Eq. ([7]) implicitly defines the largest time 
scale potentially available for computing. 
Upper-bound on the computational time- One way to 
quantify the loss of quantum information to the envi- 
ronment is through the trace distance 2] between the re- 
duced density matrix p f> (T) and the ideal density matrix 
po: D(pr{T), po) = ifr \pr{T) - po\. The trace distance 
indicates how hard it is to distinguish two density matri- 
ces by performing measurements; hence, it is a natural 
way to quantify how well QEC protects information. Let 
us suppose that there is a criterion £)(pfl(T), po) < Dent 
for a successful computation. Our goal, then, is to eval- 
uate the time T available for computation. 

Since we expect that D(^p]^{T), poj is small, it is nat- 
ural to formulate the problem in powers of the effective 
couplings A* . For an upper bound on T, we can stop the 
perturbative expansion in second order. Though it is dif- 
ficult to evaluate D(^pfi.{T), po) in general, we can make 
some progress by considering two distinct cases. First, 
we look at an isolated logical qubit, namely, S — !■ oo. 
Second, we use the Hilbert-Schmidt norm to bound the 
trace distance and define an upper bound on T in general. 
Information lost by a single logical qubit- For a sin- 
gle logical qubit, the trace distance can be expressed 
in terms of the expe ctation values of the logical 
qubit D{pRiT),po) = ^\5a+iT)\^ + [Sa-{T)]y4, where 
Sa°'{T) = {a"{T)) — (a") and, for convenience, we 
dropped the space label. Since the largest coupling con- 
stant is in the z direction, we employ a rotation to take 
it into account nonperturbatively. First, we define the 
operator : i^^ ((n + 1)A) : - :i^^(nA):= A* A :/^(nA): 
and then rotate the evolution operator at each n'^ QEC 
period using the unitary transformation e~*'^ (tiA):(t ^ 
This rotation cancels the z component of Hj at the ex- 
pense of dressing the transverse coupling. The rotated 



interacting Hamiltonian at a time t = nA can be written 
as H}°\t) = a; Ea={±} ■■ rit) ■■ exp [-2ia : F^{t) :] a". 

We can now calculate the expectation values 6a°'{T) 
in perturbation theory in X*. This is a simple but te- 
dious calculation which we omit here. For our purposes, 
the main feature appears already at zeroth order (de- 
phasing only) . In this case, it is straightforward to show 
that {a'{T)) = (ct") and {a+{T)) = e'^'^'^^'' {a+), where 

7.(T) = (2Vi)^(A*)2^k^o^[l-cos(c.,,fer)] is the 

well-known decoherence function 141 . We thus obtain 



D{pn{T),po) = \{a+)\ [l - e-^''^ (T)] 



(8) 



By defining — 2(2^ — s^) — D, we can distinguish the 
following decoherence regimes in the long-time limit: 



,,rAN , (A*/a;o) InM, 
^^(^^^)°^^(A;M)^(.oA)W^^mW^ 

(A*A)2 {koL/2n)^'-''' 



Cz<o, 

, 0<Cz < 2z^, 
M2, C, > 2z„ 

(9) 

where M = T/A is the number of QEC steps that are 
performed. These regimes are straightforward general- 
izations of the super-Ohmic {Cz < 0), Ohmic (Cz = 0), 
and sub-Ohmic (Cz > 0) regimes. Notice that for Cz < 0, 
the trace distance will converge to a finite value Dsat- 
Equation is an exact result but we expect I?crit to be 
small. Hence, we can expand the exponential and find 
the maximum time for computation with isolated logical 
qubits. Assuming I?ciit > ^sat, we find 



M, 




(u;o/A*)2^^/CV(^oA), 



Cz = o, 

< Cz < 2zz 
S2n/koLf^-^'^^ y/D—t/{XlA), Cz > 22z, 

.(10) 

where cd,z is a dimensionless pref actor of order unit. 
Upper bound for multiple logical qubits- To find an upper 
bound on the trace distance when logical qubits are not 
isolated, we use the sub-additivity property of the square 
root function and an inequality proved in Ref. (isj , 

DHsipRiT),po) < D{pr{T),po) < 2^Dhs{pr{T),p^), 

where D h s [p b{T) , po) = jMpr{T) - Po?Y'^ is the 
Hilbert-Schmidt norm and N is the number of logical 
qubits. Following a similar procedure to that used for 
the trace distance, we can expand Dhs{pr{T), po) to 
second order in A* , 



Dhs(pr{T),Pq) 



oc 



^1 I. 



-ik-(x-y) 



(T) , (12) 



(13) 
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There are two types of W^y{T): (i) the diagonal self- 
interaction terms, and (ii) the correlation terms in which 
pairs of logical qubits interact. Both types lead to the 
same functional dependence in the sum: 



Nu^^ InM, 



Ca <0, 

(14) 

where the proportionality constant is of order 1. How- 
ever, the two types of terms lead to different onset crite- 
ria. For the self- interacting part, the different regimes 
are delineated using C^a = '^{za — Sa) — D, while for 
the correlation part, the spatial sum leads to Co — 
2{za — Sa) + Dx — D with Dx being the dimension of 
the qubit array. Note that some physical arrangements 
of qubits are more favorable than others; for instance, 
a linear architecture is more favorable than a square or 
cubic one. 

First, note that in order to apply QEC we assumed 
A*^7V < 1 [Eqs. dm) and dH])]. Second, for a given 



critical distance Dent and using Eq. (1141) , we arrive at an 
upper bound on the time available to compute due to 
each component of the environment: 



oo, 
exp 



(c^oA)-i 



_Dc: 



Ca <0, 

Ca -0, 

< Ca < Za, 

(15) 

where bD,a is a dimensionless constant. If A^; ^ ^ A, 
this result is simply related to the code distance or the 
level of concatenation: A* ~ A'^'' . 



Af(A*A) ' 



Conclusions- For how long is it possible to quantum com- 
pute? (1) The trace distance calculations give us a rule 
of thumb: For a finite computation time, the residual 
decoherence of a logical qubit after the first QEC step 
times the number of logical qubits must be a small num- 
ber, (A*)^Af <C 1. (The case of no spatial correlation 
among the physical qubits at short times. A* = 0, was 
discussed in Ref. [8].) This condition must be a factor in 
the choice of the distance of the code or concatenation 
level. (2) While the argument presented here does not 
directly address the threshold theorem, the upper bound 
on the available computational time shows that there are 
certain limits to the power of QEC. The three regimes 
that we find nicely fit the qualitative interpretation of 
resilience as a dynamical quantum phase transition [sj. 

(2.1) For (x,z < (above the "upper critical dimension"), 
the usual enunciation of the threshold theorem j(| can be 
used, and therefore it is possible to compute indefinitely. 

(2.2) For (^x,y > Zx^y (below the "lower critical dimen- 
sion" ) , correlations are so strong that the available com- 
putational time is formally zero (since it depends on the 
size of the bath, L). It is, however, conceivable that its 
strong infrared divergence may be handled by combin- 
ing dynamical decoupling and QEC methods Q. (2.3) 
Finally, between these two regions, there is a maximum 
time available to compute. This constraint must also be 
a factor in the choice of the distance of the code or con- 
catenation level. Even though the regimes fit into the 
general discussion of Ref. [8] , the definition of the upper 
critical dimension given here is not the same. The reason 
is that, we have now shown that it is possible to explicitly 
treat a dense set of qubits. 

We thank D. Lidar for useful discussions and corre- 
spondence. This work was partially supported by the 
Office of Naval Research and CNPq-Brazil. 
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